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ABSTRACT 

The  objective  of  this  study  is  to  introduce  two  estimators  of  parameter  of  Weibull  distribution,  which  depends  on 
modifying  Bayes  estimator  and  proposed  Bayes  estimator,  depending  on  the  extension  loss  function.  Then,  using  a 
simulation  study  (MATLAB  program),  to  find  the  best  one  based  on  MPE  and  MSE. 
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INTRODUCTION 

Elli  and  Rao  (1986),  estimated  the  shape  and  scale  parameters  of  the  Weibull  distribution,  by  assuming  a  weighted 
squared  error  loss  function.  They  minimized  the  corresponding  expected  loss,  with  respect  to  a  given  posterior  distribution. 
Sinha  &  Sloan  (1988),  obtained  Bayes  estimator  of  three  parameters,  Weibull  distribution  compared  to  the  posterior 
standard  deviation  estimates,  with  the  corresponding  asymptotic  standard  deviation  of  their  maximum  likelihood 
counterparts  and  numerical  examples  are  given.  In  2002,  Klaus  Felsenstein  developed  Bayesian  procedures,  for  vague  data. 
These  data  were  assumed  to  be  vague  in  the  sense  that,  the  likelihood  is  a  mixture  of  the  model  distribution  with  error 
distribution.  An  extension  of  Jeffery  prior  information  with  square  error  and  loss  function  in  exponential  distribution  was 
studied  by  Al-Kutubi  (2005). 

In  this  paper,  Alkutubi  (2009)  proposed  an  extension  of  Jeffery  prior  information,  with  a  new  loss  function  and 
then  compared  it  with  standard  Bayes,  for  exponential  distribution  to  find  the  best.  In  this  paper,  we  will  introduce  two 
estimators  to  the  parameters  Weibull  distribution,  using  modify  Bayes  estimation  and  extension  loss  function,  We  will  then 
compare  between  them,  by  a  simulation  study  to  find  the  best  one  based  on  MPE  and  MSE. 

Materials  and  Methods 


Let  tj,  t2, ... ,  tn  be  the  life  time  of  a  random  sample  of  size  fl  with  distribution  function  and  probability  density 
function.In  the  Weibull  case,  we  assumed  that  the  probability  density  function  of  the  lifetime  is  given  by 


1  exp 


The  lifetime  probability  density  function  /  (t,  6,  c)  are  regarded  as  a  conditional  probability  density  function 
f  (t\9,c)  where  the  marginal  probability  density  function  of  0  is  given  by  ,  the  extension  of  Jeffery  prior 

information. 
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d  In  /(t,6,c) 
99 2 


We  find  extension  Jeffery  prior  by  taking^ (0)  oc  [7(0)]Cl,  cq  £  7?+ 
Where  fisher  information  7  (0)  is  given  by 


(«^2)=E(M)=/:M|(,r,e_Sy 


C-l  _r2 

dt=W 


Then  7(0)  =  — ,  so  the  extension  of  Jeffery  prior  information  is  given  by 

m  =  k^ 

The  joint  probability  density  function  is: 

H(tlt ... ,  tn,  0)  =  nF=i/CtI  |6')^(6') 

knCl 


-exp 


(-m 


QTl+2Ci 

The  marginal  probability  density  function  of  (tj, ...  tn,  0,  c)  is  given  by: 

fcnci(n+2c1-2)! 


P(t i,  ...,tn)  = 


n  >.n+2ci-l 


(sf=1h) 


Then  posterior  distribution  is, 

e-n~2cexV(-^j 
-  ffiti)1-’I-2c(n+2c— 2)! 


By  using  the  squared  error,  loss  function  /  —  j  ’  we  can  °^ta^n  Risk  function,  such  that 

7 7(0,0)  =£'L(0,0)  =  Jo°°c(0  -0)27r(0|t1,...tn)d0 

,  9r(9,9)  -  ,  ,  _  .  .  £t; 

Let  — =  0,  then  the  Bayes  estimator  is0i  = - 

99  J  1  n+2c1_z 

By  using  the  new  extension  loss  function  [1],  we  can  get  a  new  estimator  of  the  parameter  Weibull  distribution, 

such  that 

P(s,e)^fome^§-6)2  n(6\t1, ...  tn)d0  Solving  equation  in  above  to  get  the  second  Bayes  estimator 


RESULTS 

In  a  simulation  study,  we  have  chosen  n=30,  60,  90  to  represent  small,  moderate  and  large  sample  size,  several 
values  of  parameter  0=0.3,  1,  1.2,  and  values  of  the  new  loss  function  c  =  0.5, 1.1.  The  number  of  replications  used  was 
R=1000.  The  simulation  program  was  written  by  using  the  Matlab  program.  After  the  parameter  was  estimated,  mean 
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square  error  (MSE)  and  mean  percentage  error  (MPE)  were  calculated  to  compare  the  methods  of  estimation,  where 


X  [e-e) 

MSE  (g)-tr  and  MPE  | 


0)  = 


1000  \di  — 


The  results  of  the  simulation  study  are  summarized  and  tabulated  in  Table  1  and  Table  2,  for  the  MSE  and  the 
MPE  of  the  two  estimators  for  all  sample  sizes  and  d  values,  respectively.  It  is  obvious  from  these  tables  that,  modify 
Bayes  estimator  with  the  Jeffery  prior  information,  t^is  the  best  estimator.  But  02  in  most  cases  has  the  largest  MSE  and 
MPE 


Table  1:  The  Ordering  of  the  Estimators  with  Respect  to  MSE 


Size 

0 

C 

Theta  Hat 

9i 

e2 

0.3 

0.5 

1.1 

0.210 

0.224 

0.244 

0.225 

30 

1 

0.5 

1.1 

0.232 

0.235 

0.232 

0.238 

1.2 

0.5 

1.1 

0.226 

0.223 

0.231 

0.235 

0.3 

0.5 

0.125 

0.131 

1.1 

0.121 

0.131 

60 

1 

0.5 

0.121 

0.120 

1.1 

0.122 

0.123 

1.2 

0.5 

0.122 

0.123 

1.1 

0.121 

0.122 

0.3 

0.5 

1.1 

0.091 

0.091 

0.091 

0.091 

90 

1 

0.5 

1.1 

0.081 

0.080 

0.085 

0.085 

1.2 

0.5 

1.1 

0.081 

0.081 

0.081 

0.084 

Table  2:  The  Ordering  of  the  Estimators  with  Respect  to  MPE 


Size 

0 

Theta  Hat 

t 

0i 

02 

0.3 

0.5 

1.1 

0.210 

0.210 

0.214 

0.214 

30 

1 

0.5 

1.1 

0.211 

0.211 

0.223 

0.223 

1.2 

0.5 

1.1 

0.212 

0.212 

0.212 

0.211 

0.3 

0.5 

0.113 

0.114 

1.1 

0.113 

0.114 

60 

1 

0.5 

0.112 

0.113 

1.1 

0.112 

0.113 

1.2 

0.5 

0.111 

0.112 

1.1 

0.111 

0.112 

0.3 

0.5 

1.1 

0.004 

0.004 

0.004 

0.004 

90 

1 

0.5 

1.1 

0.004 

0.004 

0.005 

0.005 

1.2 

0.5 

1.1 

0.005 

0.005 

0.006 

0.006 
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CONCLUSIONS 

The  modify  Bayes  estimator  with  Jeffery  prior  information  0j  is  the  best  estimator,  when  compared  it  with  the 
other  estimators.  We  can  easily  conclude  that,  MSE  and  MPE  of  Bayes  estimators  decrease  with  an  increased  of  sample 
size. 
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